Diagonal hyperbolic systems with 

large and monotone data 
Part I: Global continuous solutions 

o 

o . 

^ ■ A. El Hajj 1 , R. Monneau 2 

April 13, 2009 

Abstract 

. In this paper, we study diagonal hyperbolic systems in one space dimension. Based on a new 

^| gradient entropy estimate, we prove the global existence of a continuous solution, for large and 

non-decreasing initial data. We remark that these results cover the case of systems which are 
hyperbolic but not strictly hyperbolic. Physically, this kind of diagonal hyperbolic systems 
appears naturally in the modelling of the dynamics of dislocation densities. 
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1 Introduction and main result 
1.1 Setting of the problem 



c3 



In this paper we are interested in continuous solutions to hyperbolic systems in dimension 
one. Our work will focus on solutions u(t, x) = (u l (t, ^))j=i,...,d, where d is an integer, of 
hyperbolic systems which are diagonal, i.e. 

dtu 1 + \u)d x u l = ® on (0,+oo)xR, for i = l,...,d, (1.1) 

with the initial data: 

u\0,x) = ui{x), xeR, for i = l,...,d. (1.2) 
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d d 

Here d t = — and d x = —. Such systems are (sometimes) called (d x d) hyperbolic 

systems. Our study of system ( 11.11) is motivated by consideration of models describing 
the dynamics of dislocation densities (see the Appendix, Section [5]), which is 



dtu { + A i:j u j d x u 

where (Ai)i,j=i,...,d is a non-negative s 
special case of system (|l.ll) . 

For real numbers a 1 < [3\ let us consider the box 

U = U* =1 [a i ,P 

u 



for i = 1, . . . ,d, 
symmetric matrix. This model can be seen 



3jS cl 



)i=X,...,c 



(1.3) 

which satisfies the following 



(HI) 



We consider a given function 
regularity assumption: 

the function G C°°{U), 

there exists M > such that for i = 1, d, 
|*(m)|<M for all ueU, 

there exists Mi > such that for i — 1, d, 
\<(v) -*(«)|< Mi for all v,ueU, 

where |w| = for w; = (w 1 , . . . , w; d ). Given any Banach space (E 

i=l,...,d 

rest of the paper we consider the norm on E 

IIHU d 

We assume, for all u G M d , that the matrix 



in the 



<d. 



||iu*||.e, for w = (w ] 

i=l,...,d 



,w d ) G E c 



(j(w))ij=i,...,d, where 
is non-negative in the positive cone, namely 

for all u G U, we have 



i 



(#2) 



&&*»>0 for every £ = (£i, £d) <= [0, +oo) d . 

i,j=l,...,d 



lJ=J.,...,d 

In (11.21) . each component u l of the initial data u = (uq, . . . ,Uq) is 
following property: 



assumed satisfy the 
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(F3) 



a* < uj < 

«o is non-decreasing, 



where L \ogL( 



k d x 4 G LlogL(R), 
is the following Zygmund space: 



for i 



,d, 



LlogL(R) = |/ G L 1 (IR) such that ^ |/|ln(e + |/|) < +00 J . 
This space is equipped by the following norm: 

||/lkio g L(K) = inf |/x > : JJjj- In fe + ^ < l} , 
This norm is due to Luxemburg (see Adams [U (13), Page 234]). 

Our purpose is to show the existence of a continuous solution u = (u 1 , . . . , u d ) such that, 
for i = 1, . . . , d, the function u l (t, •) satisfies (H3) for all time. 

1.2 Main result 

fu 2 \ 

It is well-known that for the classical scalar Burgers equation d t u + d x ( — 1 =0, the so- 
lution stays continuous when the initial data is Lipschitz-continuous and non-decreasing. 
We want somehow to generalize this result to the case of diagonal hyperbolic systems. In 
particular, we say that a function Uq = (uq, . . . , Uq) is non-decreasing if each component 
u l is non-decreasing for i = 1, . . . , d. 

Theorem 1.1 (Global existence of a non- decreasing solution) 

Assume {HI), (H2) and (H3). Then, there exists a function u which satisfies for all 
T > 0: 

i) Existence of a weak solution: 

The function u is solution of ( tJ.Jjj -( 0T^j . where 

11 e [L°°((0,+oo) x R)] d H [C([0,+oo);LlogL(M))] d and d x u G [L°°((0, T); L \ogL(R))] d , 

such that for a.e. t G [0, T) the function u(t, ■) is non- decreasing in x and satisfies the 
following L°° estimate: 



< \\u. 



Q\\L° 



for 



l,...,d, 



(1.4) 
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and the gradient entropy estimate: 



/ f (d x u i (t,x))dx+ [ [ V i j (u)d x u i (s,x)d x u j (s,x) dx ds <C 1} (1.5) 
J ^i=i,...,d Jo jR i,j=i,...,d 

where 

and Ci (T, d, M 1? ||wo||[L°°(R)] d ) H^a^oH [Li og L(iR)] d ) ■ 
ii) Continuity of the solution: 

The solutionu constructed in (i) belongs to [C([0, +00) x R)] d and there exists a modulus 
of continuity uj(5, h), such that for all S,h > and all (t, x) G (0, T — 5) x K, we /iawe: 



+ 5, x + — u(t, x)\ < C 2 uj(S, h) with u(5,h) = ——. — + ——, — , (1.7) 

ln(i + l) ln(i + l) 

where C 2 (T, d, M ,Mi, ||mq||[l<»(r)]<*, \\d x u \\ [LlogLmd ) . 

The key point to establish Theorem 11.11 is the gradient entropy estimate (11.51) . We first 
consider the parabolic regularization of the system (11.11) and we show that the smooth 
solution admits the L°° bound (11.41) and the fundamental gradient entropy inequality 
(11.51) . Then, these a priori estimates will allow us to pass to the limit when the regu- 
larization vanishes, which will provide the existence of a solution. Let us mention that 
a similar gradient entropy inequality was introduced in Cannone et al. |5| to prove the 
existence of a solution of a two-dimensional system of two coupled transport equations. 

Remark 1.2 Remark that assumption (H2) implies that the second term on the left 
hand side of 111.5]) is non-negative. This will imply the LlogL bound on the gradient of 
the solutions. 

Up to our knowledge, the result stated in Theorem 11.11 seems new. In relation with our 
result, we can cite the paper of Poupaud [24], where a result of existence and uniqueness 
of Lipschitz solutions is proven for a particular quasi-linear hyperbolic system. 

Hyperbolic systems ( 11.11) in the case d = 2 are called strictly hyperbolic if and only if 
we have: 

\u\u 2 ) <\u\u 2 ). (1.8) 

In this case, a result of Lax [T9j implies the existence of Lipschitz monotone solutions of 
fll.ip - fll.2p . This result was also extended by Serre (25j Vol II] in the case of (d x d) rich 
hyperbolic systems (see also Subsection 11.41 for more related references). Their results 
are limited to the case of strictly hyperbolic systems. On the contrary, in Theorem 11.11 
we do not assume that the hyperbolic system is strictly hyperbolic. See the following 
remark for a quite detailed example. 



4 



Remark 1.3 (Crossing eigenvalues) 

Condition (1.8\) on the eigenvalues is not required in our framework ( Theorem ! Here 
is a simple example of a (2 x 2) hyperbolic but not strictly hyperbolic system. We consider 
solution u = (u 1 ,^) of 

{d t u l + cos^^dxU 1 = 0, 
d t u 2 + u l sin(u 2 )d x u 2 = 0, 

Assume: 

i) oo) = 1, u 1 (+oo) = 2 and d x v} > 0, 

ii) u 2 (-oo) = —\, u 2 (+oo) = | and d x u 2 > 0. 

Here the eigenvalues A 1 (m 1 ,m 2 ) = cos(u 2 ) and A 2 (-u 1 ,m 2 ) = u l sin(u 2 ) cross each other at 
the initial time (and indeed for any time). Nevertheless, we can compute 

ii ( l 2\\ ( —sin(u 2 ) \ 

{j[U ,U )) id= l, 2 - ^ m 1 cos(m2) J , 

which satisfies (H2) (under assumptions (i) and (ii)). Therefore Theorem \l.l\ gives the 
existence of a solution to $1.S\) with in particular (i) and (ii). 

Remark 1.4 (A generalization of Theorem ll.ll ) 

In Theorem \l.l\ we have considered a particular system in order to simplify the presen- 
tation. Our approach can be easily extended to the following generalized system: 

d t v} + i {u,x,t)d x u i = h i {u,x,t) on (0,+oo)xR, for i = l,...,d, (1.10) 

with the following conditions: 

- 1 G W 1,co {U x K x [0, +oo)) and the matrix ( l j(u, x, jd is positive in the positive 

cone for all (u,x,t) £ U xlx [0, +oo) (i.e. a condition analogous to (H2)). 

- ti e W l ^(U x R x [0, +oo)), d x h} > and > for all j ^ i. 

Let us remark that our system fll.ll) - (|1.2p does not create shocks because the solution 
(given in Theorem 11.11) is continuous. In this situation, it seems very natural to expect 
the uniqueness of the solution. Indeed the notion of entropy solution (in particular de- 
signed to deal with the discontinuities of weak solutions) does not seem so helpful in 
this context. Even for such a simple system, we then ask the following: 

Open question: Is there uniqueness of the continuous solution given in 
Theorem CP] ? 



on (0,+oo)xM. (1.9) 
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In a companion paper (El Hajj, Monneau [11]), we will provide some partial answers to 
this question. 

1.3 Application to diagonalizable systems 

Let us first consider a smooth function u = (u 1 , . . . ,u d ), solution of the following non- 
conservative hyperbolic system: 

d t u{t, x) + F(u)d x u(t, x) = 0, u G U, iGK, t G (0, +oo), 

(1.11) 

u(x, 0) = u (x) iGl, 

where the space of states U is now an open subset of M d , and for each u, F(u) is a (dx d)~ 
matrix and the map F is of class C 1 ([7). The system (11.111) is said (d x d) hyperbolic, if 
F(u) has d real eigenvalues and is diagonalizable for any given u on the domain under 
consideration. By definition, such a system is said to be diagonalizable, if there exists a 
smooth transformation w = (w 1 (m), . . . ,w d (u)) with non-vanishing Jacobian such that 
(11.111) can be equivalently rewritten (for smooth solutions) as the following system 

dtw* + X i (w)d x w i = for i = l,...,d, 

where * are smooth functions of w. Such functions w l are called strict i-Riemann invari- 
ant. 

Our approach can give continuous solutions to the diagonalized system, which provided 
continuous solution to the original system (11.111) . 

1.4 A brief review of some related literature 

For a scalar conservation law, which corresponds to system ( 11.111) in the case d = 1 
where F{u) = h'{u) is the derivative of some flux function h, the global existence and 
uniqueness of BV solutions has been established by Oleinik [23] in one space dimension. 
The famous paper of Kruzhkov [18] covers the more general class of L°° solutions, in 
several space dimensions. For an alternative approach based on the notion of entropy 
process solutions, see for instance Eymard et al. [12] ■ For a different approach based on 
a kinetic formulation, see also Lions et al. [22] . 

We now recall some well-known results for a class of (2 x 2) strictly hyperbolic sys- 
tems in one space dimension. This means that F{u) has two real, distinct eigenvalues 
satisfying (11.81 ). As mentioned above, Lax [19] proved the existence and uniqueness of 
non-decreasing and smooth solutions for diagonalized (2 x 2) strictly hyperbolic systems. 
In the case of some (2 x 2) strictly hyperbolic systems, DiPerna [HI [7] showed the global 
existence of a L°° solution. The proof of DiPerna relies on a compensated compactness 
argument, based on the representation of the weak limit in terms of Young measures, 
which must reduce to a Dirac mass due to the presence of a large family of entropies. 
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This result is based on an the idea of Tartar [27J. 

For general (d x d) strictly hyperbolic systems; i.e. where F(u) has d real, distinct 
eigenvalues 

A 1 («)<•••< X d (u), (1.12) 

Bianchini and Bressan proved in a very complete paper [3], a striking global existence 
and uniqueness result of solutions to system (11.111) . assuming that the initial data has 
small total variation. This approach is mainly based on a careful analysis of the vanish- 
ing viscosity approximation. An existence result has first been proved by Glimm [15] in 
the special case of conservative equations, i.e. F(u) = Dh{u) is the Jacobian of some flux 
function h. Let us mention that an existence result has been also obtained by LeFloch 
and Liu [201 EE] in the non-conservative case. 

We can also mention that, our system (11.11) is related to other similar models in dimension 
N > 1, such as scalar transport equations based on vector fields with low regularity. Such 
equations were for instance studied by Diperna and Lions in [8]. They have proved the ex- 
istence (and uniqueness) of a solution (in the renormalized sense), for given vector fields 
in ^((O, +oo); W^(R N )) whose divergence is in ^((O, +oo); L°°(R N )). This study was 
generalized by Ambrosio [2], who considered vector fields in -^ 1 ((0, +oo); BVi oc {R N )) 
with bounded divergence. In the present paper, we work in dimension N — 1 and prove 
the existence (and some uniqueness results) of solutions of the system (ll.ip - fll.2l) with 
a velocity vector field l (u), i = 1, . . . , d. Here, in Theorem 11.11 the divergence of our 
vector field is only in L°°((0, +oo), L logL(IR)). In this case we proved the existence 
result thanks to the gradient entropy estimate (11.51) . which gives a better estimate on 
the solution. 

Let us also mention that for hyperbolic and symmetric systems in dimension N > 1, 
Garding has proved in [13] a local existence and uniqueness result in C([0, T); H "($&. )) PI 
C^Qt^T);^- 1 ^)), with s > f + 1 (see also Serre [25l Vol I, Th 3.6.1]), this result 
being only local in time, even in dimension N — 1. 

1.5 Organization of the paper 

This paper is organized as follows: in Section El we approximate the system (11.11) . by 
adding the viscosity term (ed xx u l ). Then we show a global in time existence for this 
approximated system. Moreover, we show that these solutions are regular and non- 
decreasing in x for all t > 0. In Section [3J we prove the gradient entropy inequality 
and some other e-uniform a priori estimates. In Section we prove the main result 
(Theorem 11.11) passing to the limit as e goes to 0. Finally, in the appendix (Section [5]) , 
we derive a model for the dynamics of dislocation densities. 
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2 Local existence of an approximated system 

The system (11. ip can be written as: 

d t u + {u)od x u = 0, (2.13) 

where u := (u) = and moj; is the "component by component product" 

of the two vectors u = (u 1 , . . . ,u d ), v = (v 1 , . . . ,v d ) G M. d . This is the vector in R d 
whose coordinates are given by {uov) 1 := u l v l . We now consider the following parabolic 
regularization of system (I2.13K . for all < e < 1: 

d t u £ + (u £ ) o d x u £ = ed xx u £ 

(2.14) 

u £ (x, 0) = Uq(x), with Uq(x) := u * r) e (x), 
d 2 

where d xx = —-^ and r\ e is a mollifier verify, T]J-) = -rii 1 ), such that r\ e C?°(R) is a 
non-negative function satisfying Lr] = 1. 

Remark 2.1 5?/ classical properties of the mollifier (r] E ) e and the fact that u £ G 
[L°°(R)] d , t/ien w e [C°°(M)] a! fl [^ 2 '°°(R)] d . Moreover using the non-negativity of rj e , 
the second equation of {2.1$ we get that 



IK' IU°°(R) < IpoIU°°(k)) / or i = l, 
and (H3) also implies that u £ Q is non- decreasing. 

The following theorem is a global existence result for the regularized system ( 12.141) . 

Theorem 2.2 (Global existence of non-decreasing smooth solutions) 

Assume (HI) and that the initial data u £ Q is non-decreasing and satisfies u £ Q G [C°°(R)] d n 
[W 2 '°°(R)] d . Then the system fitt® , admits a solution u £ G [C°°([0, +oo) x R)] d D 



[iy 2 '°°((0,+oo) x R)] d such that the function u £ (t,-) is non- decreasing for all t > 0. 
Moreover, for all t > 0, we have the a priori bounds: 

ll« e, *(*» -)IU°°(R) < ||«o' t |U'»(R)» f or i = l,---,d, (2.15) 



ll 9 * uS '1L<»([0,+oo);Ll(R)) " 2 H M i |U-W, « = 1, - . - , d. (2.16) 

The lines of the proof of this theorem are very standard (see for instance Cannone et 
al. [5] for a similar problem). For this reason, we skip the details of the proof. First of 
all we remark that the estimate ( 12 . 1 5H is a direct application of the Maximum Principle 
Theorem for parabolic equations (see Gilbarg-Trudinger [HI Th.3.1]). The regularity 
of the solution follows from a bootstrap argument. The monotonicity of the solution 
is a consequence of the maximum principle for scalar parabolic equations applied to 
w £ = d x u £ satisfying 
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d t w £ + (u £ ) o d x w £ + d x ({u £ )) o w e = ed xx w £ . 
Since d x u £ > this implies easily the second estimate (12.161) . 

3 ^-uniform a priori estimates 

In this section, we show some e-uniform estimates on the solutions of system (12.141) . 

Before going into the proof of the gradient entropy inequality defined in (11.51) . we an- 
nounce the main idea to establish this estimate. Now, let us set for w > the entropy 
function 

f(w) = w lnw. 

For any non-negative test function (p G C^QO, +oo) x E), let us define the following 
"gradient entropy" with w l := d x u % : 



S(t)= [ yKV) ( E />'(*'•))) dx 
\i=l,...,d J 



It is very natural to introduce such quantity S(t) which in the case ip = 1, appears 
to be nothing else than the total entropy of the system of d type of particles of non- 
negative densities w t > 0. Then after two integration by parts, it is formally possible 
to deduce from (11. ID the equality in the following gradient entropy inequality for all t > 



dS(t) 
dt 

with the rest 



Jr \ij=i,...,d J 



R(t) = J I {dtp) ( E /wl + (3*0 f E 7M j I 



where we do not show the dependence on t in the integrals. We remark in particular 
that this rest is formally equal to zero if <p = 1. 

To guarantee the existence of continuous solutions when e = 0, we will assume later 
(H2) which guarantees the non-negativity on the second term of the left hand side of 
inequality (13. 1TD . 

Coming back to a rigorous statement, we will prove the following result. 
Proposition 3.1 (Gradient entropy inequality) 

Assume (HI) and consider a function u 6 [L°°(R)] d satisfying (H3). For any 
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< e < 1, we consider the solution u e of the system given in Theorem 

\2.2\ with initial data u e = u * rj £ . Then for any T > 0, there exists a constant 
C (T, d, Mi, ||«o||[L-»(R)] d 7 ||^Mo||[Lio g L(iR)] d ) such that 



S(t) 



[ [ ^{u^w^w^ <C, with S{t)= [ f(w £ '%t r ))dx. (3.18) 

J ° jR i,j=l,...,d i=l,...,d 



where f is defined in 111. 6]) and w £ = (iu e,J )j=i,...,d = d x u £ . 

For the proof of Proposition 13.11 we need the following technical lemma: 

Lemma 3.2 (LlogL estimate) 

Let (r7 e ) eG ( 0i i] be a non-negative mollifier satisfying J R rj e = 1, let f be the function de- 
fined in jil.bV and h G L 1 (IR) be a non-negative function. Then 



i) J f(h) < +oo if and only ifh G LlogL(R). Moreover we have the following estimates: 

/ f(h) < 1 + |HUiogL(R) + II^IUhk) In (1 + ||/i|Uio S i(K)) , (3.19) 
Jr 

|Lio g L(M)<l + / /(/i) + ln(l + e 2 )||/i|| L i (R) . (3.20) 
Jr 



ii) Ifhe LlogL(R). then for every e G (0, 1] the function h e = h * r] e G LlogL(R) and 
satisfies 

H^ellxiogL(K) < C||^IUiogL(R) an d \\ h — h e \\ L log i(R ) — > as e — > 0, 
where C is a universal constant. 
Proof of Lemma 13.2b 

The proof of (i) is trivial. To prove estimate (13.191) . we first remark that, for all h > 
and n G (0, 1], we have 

h\n(h) + - J l {h >i } < h\n(h + e) < hln(e + fih) + | ln(//)|/i. 

We apply this inequality with u = — r— and integrate, we get 

maxfl, \\h Wdoklcs)) 



[ f(h) < - f fih\n(e + fih) + \\n^)\\\h\\ LH 
Jr t 1 Jr 



< - + | ln^m/illiim), 
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where we have used the definition of ||/i||,LiogZ,(R)- This gives (13. 191) using the fact that 
1 

/" > [7T7| ■ 

1 + ||"|UlogL(R) 

To prove (13.201) . we remark that, for h > \, we have e < e 2 h and 

hln(e + h)< h\n{h) + Mn(l + e 2 ) < f(h) + Mn(l + e 2 ). 
However, for < fa < -, we have in particular 

h\n(e + h) < hln(l + e 2 ). 



h\n(e + h) < / /(/ i ) + ln(l + e 2 )||/ i || i i ( 



Therefore 



From the definition of ||/i||ii gi(R), we deduce in particular (13.201) . For the proof of (ii) 
see Adams [H Th 8.20]. 

□ 

Proof of Proposition 13.1b 

First we want to check that S(t) is well defined. To this end, we remark that if w > 0, 
then 

0<f(w)<-l {w>1} +w\n(l + w). 

Which gives that 

[ f(w) < IMUw ln i 1 + IMU°°(R)) + / _1 {«;>i} < IMUw ( l + In (l + ||HU°°(R))) • 
Jr Jr e 

Now by Theorem E21 we have d x u £ = w £ E [L°°((0, +oo); L 1 (R))] d n [VF 2 '°°((0, +oo) x 
R)} d . This implies that S G L°°(0, +oo). We compute 



is® = I E /V'*)(^ M ), 

« .•_ 1 J 



Jl J 2 



/ \ / 

J R i=i,...,d J R i=i,...,d 



"<w e > 1 ) 
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Remark that these computations (and the integration by parts) are justified because 
on the one hand w e, \ its derivatives and 1 are bounded, and on the other hand w 6 ' 1 is 
in L°°((0, +00); L 1 (IR)). We know that J2 < because / is convex. To control J\, we 
rewrite it under the following form 



Ji = [ E 'WK'OQ 
• / *i=i.....d 



where 



Then, we deduce that 



x — - if x > 1/e, 
if 0<x<l/e. 



jR i=i,...,d 
= -f E M«> £ '^(^), 



Jll ^12 



^ M i,j=l,...,d ^ W i,j=l,...,d 

We use the fact that — x| < ^ for all x > and (HI), to deduce that 

|Jl 2 | < -dM 1 || [ J L°=((0,+oo),i 1 (R))] d 

2 



< -dMi||iio||[.L«>(R)]d :— Co(|po||[L oo (H)] <i ) d, Mi) 



where we have use 



lh £ '1L-((o >+ oc),Li( R) ) ^ 2\K\\L°°m, for z = l,...,d, (3.21) 
which follows from Remark 12.11 and Theorem I2.2L Finally, we deduce that 

J f S(t) <Jll + Jl2 + J 2 

< J11 + C . 

Integrating in time on (0, t), for < t < T, we get that, there exists a positive con- 
stant C (T,d,Mi, ||wo||[L°°(]R)] d , ||<9xWo||[LiogL(R)] d ) which is independent of e by (13 . 1 9H and 
Lemma [3721 (ii) such that 



S(t) + 



ft J2 W)w £ ' j w e ^ < C T + 5(0) < a 
Jo Jr ■ ,_, 1 ' 



i,j=l,...,d 
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□ 



Lemma 3.3 {W~ 1,1 estimate on the time derivative of the solutions) 

Assume (HI) and that the function uq E [L°°(R)} d satisfies (H3). Then for any < e < 
1, the solution u e of the system j2.14\ ) given in Theorem \2.2\ with initial data Uq = uo*r] e , 
satisfies the following e-uniform estimate for all T > 0; 

||^t'U £ ||[ L 2(( 0)T ) ;M /-l,l(K))] d < C|| M 0||[Loo( R )]d. 

where C = C(T,M ) > and W'^R) is the dual of the space W 1 '°°(R). 
Proof of Lemma 13.3b 

The idea to bound d t u £ is simply to use the available bounds on the right hand side of 
the equation (12 . 141) . We will give a proof by duality. We multiply the equation (12.141) 
by <p E [L 2 ((0,T),iy 1 ' oo (R))] d and integrate on (0, T) x R, which gives 



d t u e 



(0,T)xI 



(0,T)xI 



h 



2 „.£ 



■d xx u 



<p ■ ((u e )od x u e 



(0,T)xI 



We integrate by parts the term Ji, and obtain: 



\h\ < 



d T (b ■ d x u e 



(0,T)xl 



< w, 



'xY\\{L 2 {{0,T),L<= 

< 2T5 



\d x u e 



(3.22) 



\lL 2 ((0,T),W 1 >™(R))] d \\ U 0\\{L°°(m.)] d i 

where we have used inequality (13.211) . Similarly, for the term I2, we have: 



0llS y ll[L2(( 0)T ) )£ , o 



\d x u £ 



l[L 2 ((0,T),L ] 



(3.23) 

— 2T 2 ik/o 1 1 M 1 1 [^oc ( M )j d 1 1 1 1 [L 2 ((0,T),W 1 > oo (M.))] d - 

Finally, collecting (13.221) and fl3.23[) . we get that there exists a constant C = C(T,M Q ) 
independent of < e < 1 such that: 



/ < 


b ■ d t u e 


J(0,T)xR 





< C\\u \ 



[L° 



[^((O.T),^ 1 



which gives the announced result. 



□ 



Corollary 3.4 (e-uniform estimates) 

Assume (HI) and that the function u E [L°°(R)] satisfies (H3). Then for any < e < 
1, the solution u e of the system \2.1J$ given in Theorem \2.2\ with initial data u £ Q = uo*r] £ , 
satisfies the following e-uniform estimate for all T > 0: 



ll^ M£ |l[L° o ((0,+oo),L 1 (M))] d + ll Me |l[L°°((0,+oo)xR)] d + ll^t^ll^^O.T)^- 1 ' 1 

where C = C(T,M , \\u \\ [L<x>{R)] d) . 



< C, (3.24) 
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This Corollary is a straightforward consequence of Remark 12.11 Theorem 12.21 estimate 
(13.211 ) and Lemma [3731 



4 Passage to the limit and proof of Theorem 11.1 



In this section, we prove that the system (ll.ll) - (11.2l) admits solutions u in the distribu- 
tional sense. They are the limits of u £ given by Theorem 12.21 when e — > 0. To do this, we 
will justify the passage to the limit as e tends to in the system (12.141) by using some 
compactness tools that are presented in a first subsection. 

4.1 Preliminary results 

First, for all open interval I of M, we denote by 

L log = {/ e L\I) such that J\f\hx{e + \f\) < +oo 



Lemma 4.1 (Simon's Lemma) 

Let X , B , Y be three Banach spaces, such that we have the following injections 
X B with compact embedding and B <^-> Y with continuous embedding. 
Let T > 0. If {u £ ) e is a sequence such that, 

||w £ ||l°°((0,T);X) + ll^i M£ |lL9((0,T) ; y) — ^ 

where q > 1 and C is a constant independent of e, then (u £ ) £ is relatively compact in 
L P ((0,T);B) for alll<p<q. 

For the proof, see Simon [26], Corollary 4, Page 85]. 

In order to show the existence of a solution to system (II. ID in Subsection 14.21 we will ap- 
ply this lemma to each scalar component of u £ in the particular case where X = W 1 ' 1 ^), 
B — L\I) and Y = W^I) := (Wo' 00 (I))'. 

We denote by K exp (I) the class of all measurable function u, defined on J, for which, 

^"(e |n| - 1) < +oo. 

The space EXP(I) = : \i > and u G K exp (I)} the linear hull of K exp (I). 
This space is supplemented with the following Luxemburg norm (see Adams [IJ (13), 
Page 234] ): 
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\\u\\exp(i) = inf | A > : y * - lj < 1 
Let us recall some useful properties of this space. 

Lemma 4.2 (Generalized Holder inequality, Adams p], 8.11, Page 234]) 

Let h E EXP(I) and g E LlogL(I). Then hg E L 1 ^), with 

\\hg\\ L i(i) < 2||/i||£XP(/)||fl , ||LiogL(/)- 
Lemma 4.3 (Continuity) 

Let T > 0. Assume that u E L°°((0, +oo) x E) swc/t tfjcrf 

||^a;'w||L oo ((0,T);LlogL(R)) + ||^w||L°°((0,T);LlogL(lR)) < C*2 

T%en iftai /or a// 5, h > anc? a// (£, x) G (0, T — 5) x R, we have: 
\u(t + 6,x + h) -u(t,x)\ < 6C 2 ( T—p^- r+ ' 



Proof of Lemma 14. 3b 

For all h > and (t, x) E (0, T) x R, we have: 



|u(t, x + h) — u(t, x)\ < 



x+h 



d x u(t,y)dy 



< 2||l||EXF(a:,a;+?i)||^w(t, •) IUlogZ(a:,a:+fc)> 
- 2 l n (I + 1) ll^w|U~((0,T);Llo g L(R)), 



(4.25) 



< 2C 2 - ' 



where we have used in the second line the generalized Holder inequality (Lemma 14.2ft . 
Now, we prove the continuity in time, for all 5 > and (t, x) E (0, T — 5) x R, we have: 

+ 8, x) — u(t, x)\ 

"X+S 

\u(t + 5, x) — u(t, x)\dy, 

Ki K 2 K 3 



, * , , A , , 



"X+S i>x+S i>x+S 

</ \u(t + 5,x) -u(t + S,y)\dy,+ I \u(t + 8,y) - u(t,y)\dy,+ \u{t,y) - u{t,x)\dy . 

X J X 
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Similarly, as in the last estimate (14.25ft . we get that: 

rx+S rx+5 

Kt + Ks <5 \d x u(t + 5,y)\dy,+5 \d x u{t,y)\dy, 

J X J X 



<4C 2 



Now, we use that d t u is bounded in L°°((0, T); L logL(IR)), to obtain that: 

"X+5 rt+8 



K, < 



X J t 



\d t u(s,y)\ds dy, 



< 2S\\1\ 



EXP(x,x+5) || Ut «|| L°° ((0,T);L log L( 



< 2G 



6 



2 U7T 



Hi + 1) 



Collecting the estimates of Ki, K 2 and K 3 , we get that: 



\u{t + S,x)-u(t,x)\ < ~{K X + K 2 + Kz) <6C 2 — j 



ln(| + l) 



This last inequality joint to (14.251) implies the result. 



□ 



4.2 Proof of Theorem [TTTi 

The authors would like to thank T. Gallouet for fruitful remarks that helped to simplify 
of the proof of Theorem 11.11 Before to prove Theorem 11.11 we first prove the following 
result. 

Theorem 4.4 (Passage to the limit) 

Assume that u £ is a solution of system j2.14\ ) given by Theorem \2.2\ with initial data 



tin 



Mo * Ve where uq satisfies (H3). If we assume that for all T > 0, there exists a 



constant C > independent on e, such that: 

\\d x u e 



< C 



(4.26) 



[L°°((0,T);LlogL( 

then up to extract a subsequence, the function u e converges, as e goes to zero, to a 
function u weakly-k in [L°°((0, +oo) x Moreover, u is a solution of (CH]j-(D3^, 

and satisfies 



\u\ 



[L°°((0,+oo)xR)] d — ll M o||[ L oo( R )]d, 
\ ll^r M ll[L°°((0,T);LlogL(R))] ti — ^' 



I ll^l 



[L°°((0,T);LlogL 



d < MnC, 
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and u(t, •) is non- decreasing in x, for all t > and satisfies 



||w*IU»((0,+oo);£i(R)) < 2||uolU«CR) f or « = 1, ■ ■ ■ , (4-27) 
Proof of Theorem 14.4b 

Step 1 (w solution of (11.11) ): First, we remark that by estimate (13.241) we know that 
for any T > 0, the solutions u e of the system ( 12.141 ) obtained with the help of Theorem 
12.21 are e-uniformly bounded in [L°°((0, +oo) x R)] d . Hence, as e goes to zero, we can 
extract a subsequence still denoted by u £ ', that converges weakly-* in [L°°((0, +oo) x 
to some limit it. Then we want to show that u is a solution of system (11. IB . Indeed, since 
the passage to the limit in the linear terms is trivial in [D'((0, +oo) x M.)] d , it suffices to 
pass to the limit in the non-linear term 

(u £ ) o d x u e . 

According to estimate ( 13.241) we know that for all open and bounded interval / of R 
there exists a constant C independent on e such that: 

W u£ \\[L^((0,T);W 1 ' 1 (I))] d + H^*' Ue |l[L 2 ((0,T);W- 1 ' 1 (/))] d — ^ 

From the compactness of W X (I) c — > ^(1), we can apply Simon's Lemma (i.e. Lemma 
1441) . with X = \W 1 ' 1 (I)] d , B = [L 1 (I)] d and Y = [W- 1 ' 1 (I)] d , which shows in particular 
that 

u £ is relatively compact in [^(((^T) x I)] d . (4.28) 
Then, we can see that (up to extract a subsequence) 

(u 6 ) — > (u) a.e. 

Moreover, from Lemma l4~2l similarly as in (14.250 . we can get, for all t G (0, T) the 
following estimates: 



d x u £ (t,y)dy 



< 2C 



ln(^ + l) : 



where C is given in (14.261) . By the previous estimate and the fact that {u e ) is uniformly 
bounded in [L°°((0, +oo) x Wj] d and converges a.e. to (it), we can apply the Dunford- 
Pettis Theorem (see Brezis pj Th IV. 29]) and prove that 

(u e ) o d x u £ — > (u) o d x u 

weakly in [L 1 ((0,T) x I)] d . Because this is true for any bounded open interval I, then 
we can pass to the limit in (12.141) and get that, 

d t u + (u) o d x u = in D'((0, +oo) x R). 
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Step 2 (i priori bounds): By weakly-* convergence and from the fact that L°°((0, T); L log L(JNL)) 
is the dual of L 1 ((0, T); i? exp (IR)) (see Adams [1] for the definition of the Banach space 
E exp (M.)), we can check that u satisfies the following estimates: 

ll^a; U ll[£°°((0,T);£log£(R))] d < ^^J"^ II II [L°°((0,T);L log L(M))] d — ^' 

H' U ll[L° o ((0,+oo)xM)] d - ^^Jq 1 ^ ll Me H[L° c ((0,+oo)xR)] d - ll M o|| [L°c(R)] d • (4.29) 

Thanks to these two estimates, we obtain that 

\\dt U \\[L°°((0,T);L\o S L(R))] d — II ( U ) ° ^ U II [L°° ((0,T);L log L(R))] d 

< M \\d x u\\ [LOC( ( ^. LlogL(K ^d < M C. 

Moreover (14.271) follows from (14.291) and the fact that u(t, •) is non-decreasing in x (as it 
was the case for u £ ). 

Step 3 (Recovering the initial data): Now we prove that the initial conditions ( 11.2ft 
coincides with u(0, ■). Indeed, by the e-uniformly estimate given in Corollary 13.41 we 
can prove easily that, we have 

\\ u6 (t) ~ u o\\[W- 1 ' 1 (R)] d <Cti. 

Then, we get 

\\u(t) — ■Uo||[V7-i,i(R)] d — \\ u — u o\\[L°°((o,t);W- 1 ' 1 (R))] d 

< liminf \\u £ - «oll [I/ o 0((0it);W ,-i,i( R ))]<* < Ct*, 

where we have used the weakly-* convergence in L°°((0, i); W /_1,1 (]R)) in the second line. 
This proves that u(0, ■) = u in [D'(IR)] d . 

□ 

Proof of Theorem II. lb 

Step 1 (Existence): Remark that by assumption (H2) and estimate (13 . 1 8j) . we deduce 
from (I3.20P joint to (I3.2ip that, the solution u e given by Corollary 13.41 satisfies the 
following estimate: 

11*9x^11 [ LO o(( 0iT ) ;LlogL ( K)) ]d < C, (4.30) 

where C = C (T,d,M 1 , ||uo||[.l°°(r)]«i, ||o> a! Wo||[iiogi,(K)]<»)- Now ? we a PPly Theorem S3 to 
prove that, up to extract a subsequence, the function u £ converges, as e goes to zero, to 
a function u weakly-* in [L°°((0, +oo) x IR)] d , with u is beeing solution to ( ll.ll) - (11.2l) . 
Moreover, from Lemma [4731 we deduce that the function u satisfies the continuity esti- 
mate (HTTP . 



18 



Step 2 (Justification of (11.51) ): Let 

r^(« £ ) = \ &(w £ ) + 3 M)), for i,j = i,...,d, 

w e = d x u e . 

For a general matrix T, where *r = T > 0, let us introduce the square root B = vT of 
T, uniquely defined by 

l B = B > and B 2 = T. 

Remark that for non-negative symmetric matrices, the map T i — > ^/T is continuous. 
Then we can rewrite 



^{u^W^W^ = 

i d • / " 

where C is given in (13.181) . Therefore 



x/rV) 
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y/T(u e )w e -> g weakly in [L 2 ((0, t) x R)f 



Applying the same argument as in Step 1, of the proof of Theorem 14.41 for the conver- 
gence of (u £ ) od x u £ , we see that 



B ir -> V r (w)^M = 5 weakly in [^((O,*) x R)] . 
Therefore, using the weakly convergence in L 2 ((0,t) x R), we get 



[ [ J2 'J") -"'^" 1 = f f I 2 <liminf / / v^ftfj 



{u e )d x u e 



< C. (4.31) 



Remark also that for w l = d x u\ we have 



SUp / f(w*) < 1 + ||^||L°°((0,T);LlogL(R)) + 11^* II L°°((0,T);L 1 (R)) l n + 11^ II L°° ((0,T);L log L(R)) ) 

o<t<T Jr 

< 1 + H' u;, 1lL°°((0,T);LlogL(R)) + 2 II u II L°° (R) l n (l + II ^ ||£<x>((o,T);£log£(R))) : = 

< liminf ^[u; 5 '*] 

< 1 + C + 2||4|Uoc (K) ln(l + C) := C, 



where in the first line we have used (13.191) . in the second line we have used ( 14.271) . 
in the third line we have used the weakly-* convergence of w e ' 1 towards w l in 
L°°((0, T); LlogL(R)) and in the fourth line, we have used (I4.30D . Putting this result 
together with (I4311 . we get with C x = C + C. 

□ 
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5 Appendix: Example of the dynamics of dislocation 
densities 



In this section, we present a model describing the dynamics of dislocation densities. We 
refer to Hirth et al. [17] for a physical presentation of dislocations which are (moving) 
defects in crystals. Even if the problem is naturally a three-dimensional problem, we 
will first assume that the geometry of the problem is invariant by translations in the 
^-direction. This reduces the problem to the study of dislocations densities defined 
on the plane (x\,X2) and moving in a given direction b belonging to the plane (xi,X2) 
(which is called the "Burger's vector"). 

In Subsection 15.11 we present the 2D-model with multi-slip directions. In the particular 
geometry where the dislocations densities only depend on the variable x = x\ + X2, 
this two-dimensional model reduces to a one-dimensional model which is presented in 
Subsection 15.21 Finally in Subsection 15.31 we explain how to recover equation (1 1.1 ft as a 
model for dislocation dynamics with 

j=i,...,d 

for some particular non-negative and symmetric matrix A. 

5.1 The 2D-model 

We now present in details the two-dimensional model. We denote by X the vector 
X = (x\, X2) G M 2 . We consider a crystal filling the whole space M 2 and its displacement 
v = (vi,V2) '■ IR 2 — > M 2 , where we have not yet introduced the time dependence. 

We introduce the total strain e(v) = (sij(v ))j j = \ 2 which is a symmetric matrix defined 
by 

The total strain can be spitted in two parts: 

£^) = 4 + 4 with e p = Yl £ °' kuk i 

k=l,...,d 

where is the elastic strain and ■ is the plastic strain. The scalar function u k is the 
plastic displacement associated to the k-ih slip system whose matrix e^ k is defined by 
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where (b k ,n k ) is a family of vectors in M 2 , such that n h is a unit vector orthogonal to 
the Burger's vector b k (see Hirth et al. [U] for the definition of the Burger's vector of a 
dislocation) 

To simplify the presentation, we assume the simplest possible periodicity property of 
the unknowns. 



Assumption (H): 

i) The function v is Z 2 -periodic with / v dX = 0. 

J(o,i) 2 

ii) For each k = 1, . . . , d, there exists L k G M 2 such that u k (X) — L k ■ X is a Z 2 -periodic. 
Hi) The integer d is even with d = 2N and we have for k = 1, . . . , N : 

L k+N = L k , n k+N = n k , b k+N = -b k , e°> k+N = -e^ k . 

iv) We denote by T k G M 2 a unit vector parallel to b k such that r k+N = r k . We require 
that L k is chosen such r k ■ L k > 0. 

Remark in particular that the plastic strain is Z 2 -periodic as a consequence of As- 
sumption (H). The stress matrix is then given by 

k,l=l,2 

where A = (Ajj^ . fe l=l 2 , are the constant elastic coefficients of the material, satisfying 
for some constant m > 0: 

i,j,k,l=l,2 »J'=1,2 

for all symmetric matrices e = (%) i7 , i-e. such that = 



Then the stress is assumed to satisfy the equation of elasticity 



On the other hand the plastic displacement u k is assumed to satisfy the following trans- 
port equation 



d t u k = c k r k .Vu k with c k = o-ijE 
This equation can be interpreted, saying that 



0,k 
ij 

»,i=i,2 
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r k .Vu k > 0, 



(5.33) 



is the density of edge dislocations associated to the Burger's vector b k moving in the 
direction r k at the velocity c k . Here c k is also called the resolved Peach-Koehler force in 
the physical literature. In particular, we see that the dislocation density 6 k satisfies the 
following conservation law 

d t e k = div{c k r k e k ). 

Finally, for k = 1, . . . , d, the functions u k and v are then assumed to depend on (t, X) 6 
(0, +oo) x M. 2 and to be solutions of the coupled system (see Yefimov [28l ch. 5.] and 
Yefimov, Van der Giessen [29]): 



E dcr. 



3=1,2 



0~! 



Ox, 



on (0,+oo)xR 2 , fori = 1,2, 



E A ^ ( £k ^ - E £ > h ) 

jU=l,2 V k=l,...,d / 



Sij{v) 



d t u k 



1 / dvi dvj 



+ 



2 \dxj dx 



on (0,+oo) x R 2 , 
on (0,+oo) x M 2 , 



for i,j = 1,2 



E ^-e?/ I rk - Vyk 011 (°> +0 °) x 

i,3'6{l,2} 



for = 1, 



(5.34) 



where Aj^;, fixed parameters previously introduced, and the unknowns of the 

system are u = (u k )k=i,...,d and the displacement v = (vi,v 2 ). Remark also that our 
equations are compatible with our periodicity assumptions (if), 



For a detailed physical presentation of a model with multi-slip directions, we refer to 
Yefimov, Van der Giessen [29] and Yefimov [28, ch. 5.] and to Groma, Balogh [H] 
for the case of a model with a single slip direction. See also Cannone et al. [5j for a 
mathematical analysis of the Groma, Balogh model. 



5.2 Derivation of the ID-model 

In this subsection we are interested in a particular geometry where the dislocation den- 
sities depend only on the variable x = x\ + x-i. This will lead to a ID-model. More 
precisely, we make the following: 



Assumption (H'): 
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i) The functions v(t, X) and u k (t, X) — L k ■ X depend only on the variable x = X\ + %2- 

ii) For k = 1, ... id, the vector r k = (rf, r k ) satisfies r k + r 2 fc > with fi h = — ^— - — £• 

Hi) For k = 1, . . . , d, the vector L k = (L k , L%) satisfies L\ = L\ = l k . 

For this particular one-dimensional geometry, we denote by an abuse of notation the 
function v = v(t,x) which is 1-periodic in x. By assumption (H'), (Hi), we can see 
(again by an abuse of notation) that u = (u k (t, x))k=\,...,d is such that for k = l,...,d, 
u k (t, x) — l k ■ x is 1-periodic in x. 

Now, we can integrate the equations of elasticity, i.e. the first equation of (15.34j) . Using 
the Z 2 -periodicity of the unknowns (see assumption (H), (i)~(ii)), and the fact that 
£ o,k+N _ _ £ o,k ^ gee assum ption (if), {in)), we can easily conclude that the strain 



e e is a linear function of (u J — M"' +Ar )j=i,...,7v and of (^j (n J — v? +N ) dx^j 



j=l,...,N 

(5.35) 



This leads to the following Lemma 
Lemma 5.1 (Stress for the ID-model) 

Under assumptions (H), (i)-(ii)-{iii) and (H f ), {i)-{iii) and ( 15. 32) . we have 

- a : e ' 1 = A ijU j + Q - j u j dx, for i = 1, . . . , N, (5.36) 

j=i,.. .4 j=i,-4 -' {) 

where for i,j = 1, . . . , N 



(5.37) 



M,j — Aj t i and A i+ Nj — —A^j — Aij + N — —A i+ Nj + ^, 

Qi,j = Qj,i an d Qi+N,j — ~Qi,j = Qi,j+N = ~Qi+N,j+N ■ 

Moreover the matrix A is non-negative. 

The proof of Lemma EH] will be given at the end of this subsection. 

Finally using Lemma 15.11 we can eliminate the stress and reduce the problem to a 
one-dimensional system of d transport equations only depending on the function u 1 , for 
i = 1, . . . ,d. Naturally, from (15 .361) and [H'), (ii) this ID-model has the following form 
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The ID-model of the dynamics of dislocation densities: 



i d t u i +( Y A ij uj + Y ®v [ uj dx ) 
\j=i,...,d j=i,...,d J ° ) 



//'()///'+( )^ A ; jit 1 } ^ <J I //■' dx ) i) x u l = 0, on (0, +00) x R, for i = 1, . . . , d, 

(5.38) 

with from (15.331) 

dxU* > for i = 1, . . . , d. (5.39) 
Now, we give the proof of Lemma 15.11 

Proof of Lemma I5.lt 

For the 2D-model, let us consider the elastic energy on the periodic cell (using the fact 
that e e is Z 2 -periodic) 



E(u,v) = ±f Yl ^m^kidX with 4 = Y 4 V - 

^(°' 1 ) 2 i,j,k,l=l,2 k=l,...,d 

By definition of and £?•, we have for k = 1, . . . , d 

Y = -K"M- ( 5 - 4 °) 

On the other hand using (if'), (i)-(m), (with a; = xi + 22) we can check that we can 
rewrite the elastic energy as 

1 f 1 

E=- Y A ijkis e tJ e e kl dx. 

J ° i,j,k,l=l,2 

Replacing e% by its expression (I5.35H . we get: 



E -- 

2 



~ f Y ^{u? - ui +N ){u l - u i+N ) dx 

+1 Y %(/< 



{u j - u j+N ) dx) / {u l - u i+N ) dx 



for some symmetric matrices Ay = Aji, Qij = Qji. In particular, joint to ( 15.401) this 
gives exactly ( 15.361) with ( 15.371) . 

Let us now consider the functions vf = u % — u l+N such that 

f w 1 dx = for i=l,. . . ,N. (5.41) 
Jo 
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From (15.320 . we deduce that 




More precisely, for alH = 1, . . . , N and for all w % G R, we set 



,i = (w* on [0, §], 
\ -w { on [i, 1], 



which satisfies fl5.41[) . Finally, we obtain that 




Because this is true for every w = (it; 1 , . . . , w N ) G M. N , we deduce that A a non-negative 
matrix. 



We refer the reader to El Hajj [9] and El Hajj, Forcadel fTUj for a study in the special 
case of a single slip direction, i.e. in the case N — 1. 

5.3 Heuristic derivation of the non-periodic model 

Starting from the model ( 15.38ll -f l539ll where for i — 1, . . . , d, the function u l {t, x) — P-x is 
1-periodic in x, we now want to rescale the unknowns to make the periodicity disappear. 
More precisely, we have the following Lemma: 

Lemma 5.2 (Non-periodic model) 

Let u be a solution of ( 15. 38\ )- f5. 39i) assuming Lemma l5J\ and u l (t, x) — l l -x is 1-periodic 
in x. Let 



□ 



u 3 s (t, x) = u J (St, Sx), 
such that, for all j = 1, . . . , d 



for a small 5 > and for j = 1, . . . 



d, 



«J(o,0 



u 3 (0, •), as 5 



0, and m j (0, ±oo) = u 



:j+N 



(0,±oo). 



(5.42) 



Then u = (M J )j=i,... )( i is formally a solution of 




on (0, +oo) x R, 



(5.43) 



where the symmetric matrix A is non-negative and d x u l > for i — 1, . . . , d. 
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We remark that the limit problem (15.430 is of type (11.1ft when = 1. In particular, there 
are no reasons to assume that this system is strictly hyperbolic in general. Neverthe- 
less, the general case fi l > can be treated with our approach developed in Theorem ll.il 

considering the entropy / £ „«/ (Wft*)) instead of / £ /(«.*(«,,)) «fa. 
Formal proof of Lemma 15.21 : 

Here, we know that u\ — 5l l ■ x is --periodic in x, and satisfies for i = 1, . . . , d 
^dtui + I A » u s + 6 % / u s dx ) d °> u s = 0> 011 (°> +°°) x K - ( 5 - 44 ) 

\i=l,...,d j=l,...,d ^° / 

To simplify, assume that the initial data Us(0, •) converge to a function u(0, ■) such that 

the function d x Us(0, •) inside the interval ( — — , — J has a support in (— i2, -R), uniformly 

\ 2o 2o J 

in <5, where R a positive constant. Because of the antisymmetry property of the matrix 
Q (see assumption (15.371) ). and because of assumption (15.421) . we expect heuristically 
that the velocity in (15.441) remains uniformly bounded as 5 — > 0. 

Therefore, using the finite propagation speed, we see that, there exists a constant C 

f-l 1 N 

independent in 5, such that d x iis(t, •) has a support on (-R — Ct, R + Ct) C ( -^r, — 
Moreover, from ( 15.421) and the fact that 

V Qij / u 3 5 dx= V Qij / (u J s - u 3 s +N ) dx, 

JO , , *r ^0 
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we deduce that 



V Qij / dx, 

3=1,-4 J " 

remains bounded uniformly in 5. Then formally the non-local term vanishes and we get 
for % = 1, . . . , d 

i 

Aj?4 + 5 Qij / i4 dx — > AyVp, as 5^0, 

j=l,..,d J=l,-,d ° J=l,...,d 

which proves that w is solution of (15.43ft . with the non-negative symmetric matrix A. □ 
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